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Background: Possible evidence for triangular flow and other "higher harmonic flows" in high 
energy nucleus-nucleus collisions has received much attention recently. It is conjectured that sev- 
eral higher harmonic flows v m may result from initial-state geometry fluctuations in A-A colli- 
sions coupled to a radially-expanding medium. But as with "elliptic flow" V2 measurements, non- 
hydrodynamic mechanisms such as jet production may contribute to other higher azimuth multipoles 
v m as biases. Purpose: Careful distinctions should be maintained between jet-related and nonjet 
(possibly hydrodynamic) contributions to v m (e.g., "nonflow" and "flow"). In this study we con- 
sider several questions: (a) To what extent do jet-like structures in two-dimensional (2D) angular 
correlations contribute to azimuth multipoles inferred from various v m methods? (b) If a multipole 
element is added to a 2D fit model is a nonzero amplitude indicative of a corresponding flow com- 
ponent? and (c) Can 2D correlations establish the necessity of nonjet contributions to some or all 
higher multipoles? Method: Model fits to 2D angular correlations are used to establish the origins 
of azimuth multipoles inferred from ID projections onto azimuth or from nongraphical numerical 
methods. Results: We find that jet-like angular correlations, and specifically a 2D peak at the an- 
gular origin consistent with jet production, constitute the dominant contribution to inferred higher 
multipoles, and the data do not require higher multipoles in isolation from the jet-like 2D peak. 
Conclusions: Inference of "higher harmonic flows" results from identifying certain nominally jet- 
like structure as flow manifestations through unjustified application of ID Fourier series analysis. 
Although the peak structure at the angular origin is strongly modified in more-central collisions 
some properties remain compatible with relevant pQCD theory expectations for jet production. 

PACS numbers: 25.75.Ag, 25.75.Bh, 25.75.Ld, 25.75.Nq 
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I. INTRODUCTION 



Nuclear collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) have been interpreted within a hydrody- 
namic (hydro) context in terms of formation of a thermal- 
ized, flowing partonic medium with small viscosity [l|, Q ■ 
But alternative analysis of spectrum and correlation data 
has revealed a minimum-bias jet contribution whose vari- 
ation with A-A centrality and collision energy seems to 
conflict with hydro expectations 0-0] • ^ n particular, a 
large-amplitude same-side (SS) 2D peak centered at the 
origin in angular correlations on relative pseudorapidity 
r\ and azimuth 0, expected as a component of jet manifes- 
tations, persists even in central Au-Au collisions, albeit 
the SS peak is elongated on rj relative to a nominally 
symmetric jet cone [H, H[ ■ 

In response the jet interpretation of the SS 2D peak has 
been challenged in the literature. One strategy was based 
on introduction of "triangular flow" to angular correla- 
tions. According to a conjecture in Ref. [8| triangular 
flow should correspond to "triangularity" of the event- 
wise initial-state A-A overlap geometry. Triangular flow 
would explain both the SS peak and AS double peak in 
more-central Au-Au collisions previously attributed 
to Mach cones [ioj . In a more- recent strategy "higher 
harmonic flows" were invoked to explain all aspects of 
the ^-elongated SS peak or "soft ridge" [11]. The vari- 
ous harmonic flows were attributed to fluctuations in the 



initial-state A-A geom etry coupled to radial expansion of 
the bulk medium (121 |13[ . In an alternative scheme the 
elongated SS 2D peak was explained in terms of glasma 
flux tubes coupled to radial flow [3, [HI ■ 

The problem of cylindrical multipoles (Fourier series 
coefficients) derived from ID projection of all 2D angular 
correlations onto ID azimuth was considered in Ref. [Tfj| . 
That study concluded that in almost all cases claimed 
higher harmonic flows (Fourier series index m > 2) are 
actually ID Fourier components of the SS 2D peak. The 
SS peak can be represented accurately by a factorized 2D 
model function 4] but not by a ID Fourier series due to 
substantial r\ variation of the peak structure, and some 
peak properties remain consistent with jet formation. 

In the present study we confront a more subtle issue: If 
a sextupole Fourier element is added to the "standard" 
2D fit model described in Ref. Q and a nonzero sex- 
tupole amplitude results, what does that imply physi- 
cally? Does a nonzero sextupole amplitude necessarily 
imply that "triangular flow" plays a role in nuclear col- 
lisions? If added multipole model elements result in sig- 
nificant changes to other fit parameters are the data sys- 
tematics reported in Ref. 0] arbitrary and misleading? 

With recently-reported 2D angular correlation data we 
demonstrate that attempts to add extra Fourier terms 
(multipole elements) to the 2D data model of Ref. [1] 
are equivalent to minor alteration of the SS 2D peak 
model, and such modifications are not systematically sig- 
nificant. Substantial variations in several parameter val- 
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ues do arise when the number of model parameters ex- 
ceeds that required by the correlation structure in the 
data. However, the parameter changes are strongly cor- 
related and are found to be equivalent to adding a narrow 
ID Gaussian to the SS 2D peak model. That addition 
does not change the combined data model significantly. 

This paper is arranged as follows: In Sec. [TT] we intro- 
duce some basics of 2D angular correlation measurement. 
In Sees. Hill and IIVI we define the standard 2D fit model 
for this study and provide some example fits to data. In 
Sees. IVl and I VII we consider ID projections of data onto 
pseudorapidity rj and azimuth 4> differences and some re- 
sults of ID Fourier analysis on azimuth. In Sec. IVIII we 
consider consequences of adding a sextupole to the stan- 
dard 2D fit model within a limited detector 77 acceptance 
and possible interpretations. In Sec. IVIIII we determine 
what azimuth multipole elements are required by 2D data 
and what are not required and may be omitted. Sees. IIXI 
and IXl present the discussion and summary. 



II. BASIC 2D ANGULAR CORRELATIONS 

This study addresses the problem of modeling two- 
particle angular correlations and physical interpretation 
of the model elements. Different modeling methods may 
suggest different physical interpretations. Some imposed 
terminology presumes specific interpretations. One can 
question whether there is a unique "best" model and 
whether physical interpretations must depend on anal- 
ysis techniques or may be uniquely determined by data. 

A. Kinematic variables 

Two-particle correlations are structures in pair- 
density distributions on six-dimensional momentum 
space (ptii Vi> 4 > iiPt2, ??2, 02)- 2D correlations on trans- 
verse momentum p t or transverse rapidity y t = hi[(pt + 
mt) / itItt] (m-ir is assumed for unidentified hadrons) are 
complementary to 4D angular correlations in the 6D pair 
momentum space, yt is preferred for visualizing correla- 
tion structure on transverse momentum. Angular corre- 
lations can be measured by integrating over the entire 
(pt,Pt) pah acceptance (minimum-bias angular correla- 
tions) or over subre gion s (e.g., "trigger-associated" di- 
hadron correlations) [lj, [l8| ■ 

Two-particle angular correlations are defined on 4D 
momentum subspace (771, 772, <f>i, $2). Within acceptance 
intervals where correlation structure is approximately in- 
variant on mean azimuth or polar angle (e.g. 77s = 771+772) 
angular correlations can be projected by averaging onto 
difference variables (e.g. 77A = 7/1 — 772) without loss of in- 
formation to form angular autocorrelations [l9l . |20| . 2D 
subspace (?7a,0a) is then visualized. The pair angular 
acceptance on azimuth can be separated into a same-side 
(SS) region (|</>a| < i"/2) and an away-side (AS) region 
(|^a| > tt/2)- The SS region includes mirajet corre- 



lations (hadron pairs within single jets), while the AS 
region includes internet correlations (hadron pairs from 
back-to-back jet pairs). 

Centrality variation of the A-A collision geometry is 
represented by the Glauber model [2l|. The model pa- 
rameters are the number of binary N-N collisions Nbin 
and the number of participant-nucleon pairs N part /2 
as functions of the fractional A-A cross section ct/o-q. 
We also define a mean participant path length v — 
2N},i n /Np ar t. The A-A impact parameter is represented 
by symbol b. 



B. Correlation measures 

Correlations can be measured with per-particle statis- 
tic Ap/^fp^J = p ((r) - 1), where Ap = p sib - p ref 
is the correlated-pair density, p si b is the sibling (same- 
event) pair density, p re f is the reference or mixed-event 
pair density, (r) is the mean sibling/mixed pair-number 
ratio, and prefactor pq — n c h/ArjA(j) ~ d 2 n c h/dr]d(j) is 
the charged-particle 2D angular density averaged over 
some angular acceptance ( A77, Acj>) 0, [2(| . Factorization 
Pref ~ Po is assumed, and symbol Ax denotes a detector 
acceptance on x. 

The per-particle quadrupole component of 2D angular 
correlations resulting from a given 1)2 analysis method 
is Aq {method} = po(b)v% {method} (b). The quadrupole 
component can also be measured with per-pair statis- 
tic Ap/p re f — > v%{method}, including total azimuth 
quadrupole f|{2} (method {2} denotes the m = 2 
quadrupole cosine component of all ID azimuth pro- 
jections). Higher multipoles are represented by tj^{2}. 
Variation of per-pair correlation measures with A-A cen- 
trality is typically dominated by a trivial l/n c h factor, 
or in the case of v m a l/y/n c h factor. 



C. Example 2D data histograms 

Figure Q] shows examples from analysis of 200 GeV 
Au-Au collisions [4]. On the left are most-peripheral 84- 
93% data that approximate N-N (p-p) collisions. On the 
right are most-central 0-5% data. These are minimum- 
bias (pt-integral) data: the only momentum cut is pt > 
0.15 GeV/c. There are no "trigger-associated" p t cuts 
applied. The basic correlation histograms have not been 
further processed (no "background" is subtracted). 

Interpretation of the N-N data structure in the left 
panel is straightforward. There is a broader same- 
side (SS) 2D peak at the angular origin attributed to 
minimum-bias jets or minijets. Superposed on that is 
a narrower 2D peak attributed to Bose-Einstein correla- 
tions (BEC) plus conversion-electron pairs. An away-side 
ID peak on azimuth (AS ridge) is attributed to back- 
to-back jet pairs. A ID peak on ?7a is attributed to 
dissociation (fragmentation) of projectile nucleons. No 
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FIG. 1: pt-integral 2D angular correlations from 200 GeV 
Au-Au collisions. Left: 84-93%-central collisions (approxi- 
mately N-N collisions). Right: 0-5%-central collisions. 

quadrupole cos(20a) component is evident in the periph- 
eral data. 

Figure [T] (right panel) includes only two of those fea- 
tures: the SS 2D peak and the AS ID peak. For those 
most-central Au-Au data the azimuth quadrupole has a 
negligible amplitude. The SS peak dominates angular 
correlations in all cases. Whereas the SS peak in N-N 
collisions is elongated 2:1 in the azimuth direction the 
SS peak in central Au-Au collisions becomes elongated 
3:1 in the pseudorapidity direction. The 77-elongated SS 
peak has been referred to as a "soft ridge." Although 
there is no formal difference between the peak shapes for 
two limiting centralities (both are 2D Gaussians) there 
is a strong preference in the heavy-ion literature to in- 
terpret the elongated SS peak in more-central collisions 
in terms of flows. In this analysis we focus on central- 
ity evolution and model descriptions of the SS 2D peak, 
and its manifestations in various "higher harmonic flow" 
measurements. 



D. Correlation structure terminology 

Measured angular correlations on (?7a,0a) include 
only a few components and can be modeled accurately by 
simple functional forms. In this study we distinguish be- 
tween components of measured angular correlations and 
elements of 2D fit models. Two of the fit-model elements 
(SS 2D peak, AS ID peak) have been interpreted in terms 
of measured and theoretical jet systematics One 
element (quadrupole) has been interpreted in terms of 
"elliptic flow," a conjectured hydro phenomenon. 

Conventional references to correlation structure may 
be confusing. In an assumed hydro context correlation 
structures are categorized as either flows or nonflows. 
But a jet interpretation for some structures is compelling. 
Jets are expected to appear in descriptions of high-energy 
nuclear collisions, at least in p-p and more-peripheral A-A 
collisions. In contrast, the systematics of nominally flow- 
related structures appear to contradict some aspects of 
hydro theory |22h24| . We therefore elect to categorize 
correlation structure as either jet-related or nonjet. We 
continue to use the same terminology for all Au-Au cen- 



tralities, although the jet interpretation for some struc- 
tures may be questioned in more-central collisions. 

Jet-related structure then denotes SS 2D and AS ID 
peaks. The former is modified in more-central Au-Au 
collisions, but the question remains whether a flow de- 
scription becomes appropriate there. Nonjet structure 
includes the BEC-electron peak, the ID Gaussian on r)A 
(projectile dissociation) and a quadrupole unrelated to 
SS 2D and AS ID peaks. Nonjet quadrupole systematics 
clearly demonstrate strong correlation with the initial- 
state A-A geometry. Evidence in data for other nonjet 
multipoles is the principal subject of the present study. 



E. Fit-model criteria for 2D histograms 

Attempts to isolate 2D or ID correlation structures 
from a (large) combinatoric background have generally 
followed one of two methods: (a) ZYAM (zero yield at 
minimum) subtraction from ID dihadron correlations on 
azimuth Q and (b) model fits to 2D histograms [!, H[ ■ 
Method (a) depends on background estimation, and in 
particular measurements of parameter V2 by nongraph- 
ical numerical methods which may include substantial 
jet-related bias [22|, |23j, |25[ . If the subtracted background 
is consequently biased the inferred jet-related structure 
may be underestimated and distorted [26j . 

Method (b) requires definition of a 2D fit model func- 
tion including several elements. A range of possibilities 
exists, some choices being motivated by physical mod- 
els or expectations for observation of certain phenom- 
ena. We prefer a mathematical model that is not moti- 
vated by physical models or expectations but does satisfy 
a few general requirements for statistical models: The 
same model should apply to all collision centralities, not 
just restricted intervals. The number of model parame- 
ters should not exceed what is required by the structure 
in the data (parsimony). Each model element should 
be established as necessary: omission of that element 
should result in a substantially worsened fit with equiv- 
alent residuals structure. A sufficient model is a sum 
of necessary elements that results in negligible residuals 
structure. A necessary and sufficient model is ideal. 



III. STANDARD 2D FIT MODEL 

All minimum-bias (pt-integral) 2D angular correlations 
from Au-Au collisions at ^/sjvjv = 62 and 200 GeV are 
observed to include three main components: (a) a same- 
side (SS) 2D peak at the origin on (j]a, 0a) well approxi- 
mated by a 2D Gaussian, (b) an away-side (AS) ID peak 
(or "ridge") on azimuth well approximated by AS az- 
imuth dipole cos(0a) and uniform on n& within a few 
percent (having negligible curvature) over the angular ac- 
ceptance of the STAR TPC (An = 2, Acf> = 2tt), and (c) 
nonjet (NJ) azimuth quadrupole cos(20a) also uniform 
on r/A to a few percent. Model elements (a) and (b) have 
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been interpreted together as representing minimum-bias 
jets or "minijets" [J, H, @]- Element (c) has been con- 
ventionally attributed to elliptic flow, a hydrodynamic 
phenomenon [27| . However, alternative interpretations 
may be possible [25j ]. 

Given that data structure the "standard" 2D model 
function for more-central Au-Au collisions is p| 0, [HJ 



bins at the origin. In more-peripheral collisions the BEC 
peak must be modeled explicitly. 



A + Ax exp ■ 



+ A D cos(</> A -7t)} +A q cos(2 a ). 



(1) 



A ID Gaussian on 7/a modeling projectile nucleon frag- 
mentation and a 2D exponential modeling Bose-Einstein 
correlations (BEC) and conversion electron pairs are 
omitted for simplicity in discussion of more-central 
Au-Au collisions where the former has negligible ampli- 
tude and the latter is very narrow, restricted to three 
histogram bins near the origin. The model fits presented 
here also include a small-amplitude (few percent) modu- 
lation of the AS dipole on i]& [e.g., panel (a) of Fig. [2]. 
Details are omitted here to simplify the presentation. 
The modulation is described in Sec. VII C of Ref. [J]. 

The definition of quadrupole measure Aq is statisti- 
cally compatible with jet- related measures A\ and Ad 
(all are per-particle measures), permitting quantitative 
comparisons between jet and nonjet quadrupole system- 
atics (e.g., comparing jet- related and nonjet contribu- 
tions to Aq). The standard fit model accurately de- 
scribes most p-p and Au-Au 2D data (with respect to 
centralities and p t cuts). In specific cases (e.g., ultra- 
central Au-Au collisions and "trigger-associated" high-p t 
cuts) certain modifications may be required to the nom- 
inally jet-related SS and AS peak models. 

For the present study the adopted model-parameter 
definitions in Eq. (JXJ) are compatible with those in Ref. [1] . 
In other, related papers (e.g., [l6|, 26]) certain parameter 
definitions (e.g., Aq and An) differ by factors 2 to be 
consistent with Fourier series definitions. 



IV. MODEL FITS TO AU-AU DATA 

As an example of model fits to 2D angular correla- 
tions we present data from 0-5% central Au-Au collisions 
reported in Ref. jjj . Similar fit quality is obtained for p t - 
integral correlations from all Au-Au collision centralities. 



A. Typical fit results, 0-5% central Au-Au 

Figure [2] illustrates typical 2D model fits to angular 
correlations. The best-fit model (a) is compared to the 
fitted data (b). In panel (c) are the fit residuals (2x 
more-sensitive scale), generally consistent with statisti- 
cal uncertainties. The BEC-electron peak has been sup- 
pressed in these data because it is confined to only three 




FIG. 2: Model fit to 0-5%-central 200 GeV Au-Au collisions: 
(a) model, (b) data, (c) fit residuals, (d) inferred jet structure. 

Panel (d) estimates the true jet contribution (SS 2D 
and AS ID peaks). For 0-5% collisions the nonjet 
quadrupole amplitude is consistent with zero (2^ . In 
more-central Au-Au collisions the AS dipole associated 
with the AS ID peak is observed to have a small modu- 
lation on ?7a, as discussed in Sec. VII-C of Ref. j4j. That 
modulation is included in the 2D fit model for this study 
and, with the fitted offset, has been subtracted from the 
data histogram in panel (b) to obtain panel (d). 



B. Evolution with Au-Au centrality and p t 

The full centrality dependence of the model elements is 
essential for proper physical interpretation of 2D angular 
correlations. Figure [3] shows the centrality dependence of 
the SS 2D peak amplitude (left panel) and AS f D dipole 
amplitude (right panel) [4[ . The dashed curves represent 
Glauber linear superposition (GLS) trends expected for 
jet production in case A-A collisions are linear superposi- 
tions of N-N collisions (i.e., transparent). The form of the 
dashed curves is determined by the Glauber model and 
amplitudes are matched to the peripheral Au-Au data. 

These nominally "jet-related" Au-Au structures are 
observed to follow the GLS reference accurately from p-p 
(N-N) up to a specific Au-Au centrality v 3 where the 
amplitude trends undergo a sharp transition (ST). The 
ST centrality corresponds to approximately 50% of the 
total cross section. Those trends and other evidence sug- 
gest that jet production is the correct interpretation of 
the SS 2D and AS ID correlation structures 0,01 ■ It is 
notable that the SS 2D peak amplitude scales with energy 
as log(y / iivw) as expected for a QCD process (factor 0.63 
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FIG. 3: Variation with Au-Au centrality for two model pa- 
rameters. Left: Same-side 2D peak amplitude. Right: Away- 
side ID peak amplitude. Glauber linear superposition (GLS) 
trends (dashed curves) correspond to jet production accord- 
ing to N-N binary- collision scaling. A sharp transition (ST) 
in these jet-related parameter trends near v — 3 is notable. 




FIG. 4: Left: Data from a p t -integral analysis of 9-18%- 
central 200 GeV Au-Au collisions from which fitted AS dipole, 
nonjet quadrupole and constant offset components have been 
subtracted, leaving the SS 2D peak. Right: Projection of the 
2D histogram in the left panel onto r]A- The dashed curve is 
a Gaussian with amplitude 0.20 and r.m.s. width 2.7. 



in this case) whereas the AS ID amplitude is unchanging 
with energy, as expected for dijet production |4|. 

If marginal pt cuts are applied (pt of one particle is 
constrained but not the other) the shape of the SS 2D 
peak remains a 2D Gaussian below p t 4 GeV/c, in 
contrast to so-called trigger-associated p t cuts for which 
both particles of a pair are constrained. In the latter 
case a large fraction of any jet may be excluded from the 
correlation structure Q. Generally speaking, evolution 
of both SS 2D and AS ID peak structures with p t cuts is 
consistent with jet expectations pjlGil- Deviations from 
the "standard" 2D fit model for extreme cases such as 
ultra-central (0-1%) collisions are discussed in Sec. IIX Dl 



V. PROJECTIONS ONTO ID SUBSPACES 

Given the simple structure of 2D angular correlations 
we focus on the properties of the SS 2D peak and its 
contribution to various applications of ID Fourier series. 



A. Isolating the SS 2D (jet-related) peak 

We have established that all 2D histograms for more- 
central Au-Au collisions (except "ultracentral" 0-1% 
data) can be expressed as the sum of an AS dipole 
(to = 1), nonjet quadrupole (to = 2) and SS 2D peak. 
The two multipoles are necessary model elements, as we 
demonstrate below. Thus, because of orthogonality of 
Fourier terms, any "higher harmonics" must be derived 
from the SS 2D peak structure. The SS 2D peak model 
may be questioned, but the data peak can be uniquely 
isolated with negligible model dependence simply by sub- 
tracting the fitted AS dipole and NJ quadrupole terms. 

Figure [4] (left panel) shows the result of subtracting 
fitted AS ID dipole and constant offset terms from the 
0-5% central Au-Au data in Fig. |2] (the narrow BEC- 
electron peak is truncated in the plot). No statistically- 



significant AS structure remains, no evidence for "higher 
harmonics" unrelated to the SS 2D peak. Similar results 
are obtained for other centralities. We conclude that the 
only source of Fourier components with to > 2 is the SS 
2D peak. 

B. Factorizing the SS 2D peak 

Because in almost all cases the SS 2D peak is the only 
significant source of higher multipoles we now consider 
the SS peak structure. We observe that for all cases the 
peak is factorizable. The SS peak can be described ac- 
curately by the product G(<Pa)H(jj^). Azimuth factor 
G{4>a) is always described by a narrow Gaussian. Fac- 
tor H(j]a) is also well approximated in most cases by a 
Gaussian whose width on t]a increases significantly with 
Au-Au centrality. In some cases (trigger-associated p t 
cuts, higher collision energies) H(r)j\) may deviate sub- 
stantially from a Gaussian function. To further pursue 
the higher-harmonics issue we construct projections of 
the SS 2D peak onto ID tja and </>a for pt-integral data. 

Figure |4] (right panel) shows a full 27r projection of the 
SS 2D peak in the left panel onto t/a- The small and nar- 
row peak is BEC plus conversion electrons. The Gaus- 
sian (dashed) curve passing through the data histogram 
is 0.20exp[-(TO\/2.7) 2 /2]. The SS peak for p t -integral 
data is well described by a single ID Gaussian on va- 
The same is true for most marginal ^-differential data. 

Figure [5] (left panel) shows the Fourier coefficients of 
a periodic array of unit-amplitude ID Gaussians as a 
function of its width [26(. The coefficients are given by 

F m (<J4,J = exp{-TO 2 a^ A /2}. (2) 

There are two relevant width regimes, for SS 2D and AS 
ID jet-related peaks, denoted by the hatched bands. In 
the former case several multipoles are relevant. In the 
latter case, especially for p t -integral data, only the m = 
1 AS dipole is relevant. If "high-p t " cuts are imposed 
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Gaussian width a+ A § 

FIG. 5: Left: Fourier coefficients for a periodic array of 
unit-amplitude ID Gaussians on azimuth as a function of 
the r.m.s. width. The hatched bands correspond to typical 
peak widths for same-side (SS) and away-side (AS) jet-related 
peaks. Right: The histogram is a projection of the SS peak in 
Fig. U (left panel) onto azimuth. The dash-dotted and dashed 
curves are the SS ID Gaussian and AS dipole from the 2D 
model fit. The bold dotted curves are the SS multipoles in- 
ferred from Eq. ([2]). 

the AS peak azimuth width may be reduced, and a ID 
Gaussian then becomes the preferred AS ID peak model. 

Figure[5] (right panel) shows the SS 2D peak from Fig.0] 
(left panel) projected onto 0a (bold solid histogram). 
The fitted SS ID Gaussian (dash-dotted curve) and AS 
dipole (dashed curve) are also shown. The Fourier com- 
ponents (azimuth multipoles) of the SS peak obtained 
from Eq. ^ are represented by bold dotted curves 
with amplitudes A x {SS}(b) = 2p (b)vl l {SS}(b) = 
Fm(^^ A )Ai iPro j(b), where Ai tPro j is the amplitude of the 
SS 2D peak when projected onto ID azimuth and letters 
X = D, Q, S, O denote cylindrical azimuth multipoles: 
dipole, quadrupole, sextupole, octupole, with pole num- 
ber 2m. The standard fit model of Eq. ([IJ provides an 
accurate description of 2D angular correlations, their ID 
projections and any multipoles derived from ID Fourier 
analysis of the SS 2D peak. 



have demonstrated that such nongraphical numerical 
methods (EP, two-particle cumulants) are equivalent to 
fitting the sum of all angular correlations projected onto 
ID c/>a with a single sinusoid cos(2</>a) [23]. Thus, such 
V2 results must include contributions from the jet-related 
SS 2D peak 0,[2^]. Several schemes have been proposed 
to distinguish "flow" from "nonflow" (i.e., some fraction 
of the jet contribution) (2?J . But alternative v-i methods 
typically rely on strong assumptions about the properties 
of statistical measures (e.g.. i>2{4|) that are actually con- 
tradicted by the data [22J, |23[ . And assumptions about 
how "nonflow" (jet structure) is distributed on r\ can also 
be questioned. The consequence may be substantial jet- 
related bias in published V2 results [H, |23[ . 



B. Dihadron correlations and ZYAM subtraction 

Dihadron correlations on ID azimuth with various 
trigger-associated pt cuts are expected to reveal jet struc- 
ture in more-central A-A collisions where event-wise jet 
reconstruction is difficult. Basic (raw) dihadron correla- 
tions include a large combinatoric background. Attempts 
to isolate jet-related correlation structure typically rely 
on subtracting an estimate of that background [30|. 

The background estimation is based on some criterion 
for determining a constant offset value and a quadrupole 
modulation amplitude derived from published vi data. 
The conventional method is denoted ZYAM (zero yield 
at minimum) subtraction. If V2 data are biased (e.g., by 
some fraction of the SS 2D peak quadrupole component), 
and/or the ZYAM offset criterion is invalid (e.g., SS and 
AS peaks overlap significantly) , the background estimate 
may be very inaccurate, leading to underestimation and 
distortion of inferred jet structure [26| . 



C. "Triangular flow" on ID <^a 



VI. IMPOSING A ID FOURIER SERIES ON A 

This section establishes context for the main subject of 
this paper presented in Sec. IVIII Fourier series analysis 
applied to projections of all 2D angular correlations onto 
ID azimuth has been coupled with interpretation of each 
sinusoid term as a flow manifestation. But Fourier series 
derived from ID projections cannot describe 2D angular 
correlations accurately, and ID projections can be rep- 
resented more efficiently with alternative peaked model 
functions. Several examples are considered below. 

A. "Elliptic flow" measurements 

The first report of V2 measurements (azimuth 
quadrupole) derived from RHIC collisions is Ref. [28[ 
where the "event-plane" (EP) method was invoked. We 



One of the consequences of ZYAM subtraction is the 
emergence of a double peak in the resulting AS correla- 
tion structure in place of the expected single broad AS 
peak (AS dipole) corresponding to back-to-back jet cor- 
relations [30]. The unexpected AS structure led to theo- 
retical conjectures about formation of Mach cones by in- 
teractions of energetic partons with a dense medium [Io| . 

More recently it was proposed that such AS structure 
may be a final-state manifestation of "triangular flow" 
arising from corresponding structure ( "triangularity" ) in 
the initial-state A-A geometry, via radial expansion of 
a bulk medium [§}. The underlying assumption of such 
triangular flow analysis is that all correlation structure 
projected onto ID azimuth represents flows described by 
the terms of a ID Fourier series representation. But any 
distribution on periodic azimuth can be described accu- 
rately by a ID Fourier series, and there is no obligation 
to interpret the Fourier terms as representing flows. 
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D. The 77 dependence of "triangular flow" 

Based on the assumption that inference of triangular 
flow from ID Fourier analysis of projected 2D angular 
correlations is a valid concept it is proposed to measure 
the 77 dependence of triangular flow by binning the differ- 
ence variable r\t\ and representing the projected azimuth 
correlations in each 77A bin by a Fourier series, with the 
rn = 3 sextupole term emphasized [Hj]. We apply that 
method to a particular case (centrality) with direct com- 
parison to standard 2D correlation analysis. In principle 
the same argument should apply to any A-A centrality. 



underlying peak is a Fourier component of the jet-related 
SS 2D data peak. The solid curve (ID Gaussian) is de- 
rived from the 2D model fit. The BEC-electron peak is 
proportionately much larger compared to the jet peak in 
this case. The BEC peak is narrower on azimuth, so the 
Fourier components for greater m are relatively larger, 
consistent with Eq. @. This exercise demonstrates that 
previously-measured 2D angular correlations, including 
jet-related structure, can predict any "triangular flow" 
result obtained as a function of 77A with a non-graphical 
numerical method. Any such Fourier terms with m > 2 
must be manifestations of the monolithic SS 2D jet peak. 




FIG. 6: Left: Data from a p t -integral analysis of 46-55%- 
central 200 GeV Au-Au collisions from which fitted AS dipole 
and nonjet quadrupole (plus constant offset) components 
have been subtracted, leaving the SS 2D peak. The narrow 
peak at the origin is Bose-Einstein correlations (BEC) plus 
conversion-electron pairs. Right: The SS azimuth sextupole 
component of the data in the left panel for each t?a bin. The 
solid curve is derived from the fitted SS 2D (jet) peak model. 

Figure [5] (left panel) shows the SS 2D peak obtained 
from 46-55%-central 200 GeV Au-Au collisions Q . The 
standard fit model consisting of AS dipole plus NJ 
quadrupole plus SS 2D peak provides an excellent data 
description. In the left panel the fitted offset, AS dipole 
and nonjet quadrupole have been subtracted from the 2D 
histogram leaving the only possible source of "higher har- 
monics" with to > 2. That centrality bin lies just below 
the "sharp transition" in angular correlations, in the in- 
terval where jet-related structure strictly follows the GLS 
scaling of N-N jet production expected for "transparent" 
A-A collisions (see Fig. [3]). The SS peak volume is pre- 
dicted by pQCD binary-collision scaling of dijet produc- 
tion Thus, a jet interpretation for the SS 2D peak is 
well supported. 

It is straightforward to convert measured SS 2D peak 
properties to v^{SS}(rjA)- First we separate the fitted 
narrow BEC-electron peak from the underlying broader 
SS 2D jet peak with their different azimuth widths. We 
convert the latter to the broader v^{SS}(t]a) distribu- 
tion on ?7a using the measured 2D peak properties and 
Eq. @, especially the fixed azimuth width. We do the 
same with the narrower fitted BEC-electron peak (with 
its much smaller azimuth width) and recombine the two. 

Figure [5] (right panel) shows the result. The broader 



E. Lower and higher "harmonic flows" 

The same procedure can be carried out with lower- 
order (to < 3) multipoles as well. Figure [7] illustrates 
the procedure applied to the SS 2D peak in Fig. [5] (left 
panel) for m = 1 (SS dipole) and 2 (SS quadrupole). The 
relative heights of jet-related and BEC-electron peaks 
depend on the different azimuth widths according to 
Eq. The right panel shows the dominant (jet- related) 
source of "nonflow" bias in conventional i>2 data. The 
hallmark of such SS 2D peak contributions is a large cur- 
vature on in contrast to the approximate uniformity 
on t]a of the nonjet (NJ) quadrupole within the STAR 
TPC acceptance. 
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FIG. 7: Left: The SS azimuth dipole component of the data 
from Fig. [6] (left panel) for each 77 a bin. Right: The SS az- 
imuth quadrupole component of the same data. The solid 
curves are derived from the fitted SS 2D peak model. 

Higher harmonic flows (to > 2) inferred from ID pro- 
jections onto azimuth have also been obtained from LHC 
data [HI, |33[ . In those studies imposed "77 exclusion cuts" 
are intended to exclude jet contributions assumed to be 
restricted to an 77 interval near the origin. But the def- 
inition of jet structure in more-central A-A collisions is 
a major question as yet unresolved. In such analysis the 
presence of the sharp transition in SS 2D peak properties 
is still evident at LHC energies. In Ref. 16| it is demon- 
strated that all trends associated with "higher harmonic 
flows" are predicted by jet-related angular correlations. 
The same multipole trends observed at 62 and 200 GeV 
are also evident at 2.76 TeV. The only difference is a com- 
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mon extrapolation factor 1.3 applied to all such trends. 

VII. A SEXTUPOLE IN THE 2D FIT MODEL 

We now confront the central question for this study. 
We have demonstrated that any "higher harmonics" in- 
ferred from ID Fourier analysis applied to projections of 
2D angular correlations onto azimuth are at least dom- 
inated by a contribution from the SS 2D peak, and SS 
peak properties are consistent with jet production in at 
least some cases. We now consider addition of a separate 
sextupole element to the standard 2D fit model moti- 
vated by a desire for increased sensitivity to "triangular 
flow." If we observe a significant nonzero sextupole am- 
plitude does that imply a nonjet sextupole component in 
the data that might indeed represent triangular flow? 

To answer that question we compare three 2D model 
configurations applied to the same data histogram: (a) 
the standard 2D model, (b) the standard fit model plus 
SS sextupole element, and (c) the standard fit model plus 
SS ID Gaussian on (f>&. The relevance of configuration 
(c) will become apparent below. With p t -integral data 
we show that if a sextupole element uniform on 77A is 
added to the standard 2D fit model nonzero sextupole 
amplitudes are indeed inferred from model fits to more- 
central Au-Au data. But other model parameters are also 
shifted substantially from their original values. The com- 
bination of changes is equivalent to a single new model 
element. The added sextupole element does not reveal 
an independent nonjet sextupole component in the data. 

A. Direct comparison of three 2D fit models 

Table Q] shows parameters from three 2D fits to the 
same data histogram from 9-18% central Au-Au colli- 
sions. Other more-central Au-Au data give similar re- 
sults. Fit parameters that do not change significantly 
between models are not relevant to this discussion and 
are omitted from the Table. The parameter labels are as 
defined in Eq. d). 

The second column (standard) shows fit results ob- 
tained with the standard 2D model. Those values are 
consistent within fit uncertainties with the results in Ta- 
ble III of Ref. jl| . The third column (std + A s ) shows fit 
results when a sextupole term ^4scos(3<^a) is added to 
the standard model and all else remains the same. The 
fourth column shows fit results when a SS ID Gaussian 
on azimuth Aid exp[— (</>a/o>a,i,d) 2 /2] is added to the 
standard fit model and all else remains the same. 



B. Model equivalence: sextupole vs ID Gaussian 

We now demonstrate that addition of a sextupole term 
As to the standard model is equivalent to modifying the 
SS 2D peak model by adding a SS ID Gaussian narrow 



TABLE I: Model parameters for three fit models: (a) standard 
model for this analysis, (b) standard model plus additional 
sextupole term As and (c) standard model plus additional 
SS ID Gaussian on <^a with amplitude Am as part of the 
SS 2D peak model. The model parameters are as defined in 
Eq. (0. Uncertainties in the second column illustrate typical 
fit uncertainties (statistical plus systematic) for each parame- 
ter. An denotes the additional (new) model parameter - As 
or Am in third and fourth columns respectively. 



parameter 


standard 


std + A s 


std + A 1D 


Ai 


0.76 ± 0.04 


0.51 


0.47 




2.3L0.3 


1.78 


1.72 


Aq 


0.18L0.01 


0.23 


0.18 


A D 


0.29L0.02 


0.175 


0.28 


An 




0.012 


0.29 



on <j>& and uniform on r]/±. The argument is based on 
results presented in Refs. [H, Ull where periodic peak 
arrays are approximated by truncated Fourier series. 

The model elements Ap (dipole) and Aq (quadrupole) 
already present in the standard 2D fit model are in effect 
parts of a (truncated) Fourier series. Adding an As (sex- 
tupole) term to the data model can extend the Fourier 
series to m = 3. The next term Ao (octupole) is typi- 
cally at the level of statistical uncertainties [1, UH ■ The 
separate sinusoids may serve both as representatives of 
distinct (nonjet) data multipoles and as parts of a Fourier 
series representing a localized (possibly jet-related) peak 
structure on azimuth [lfjl [26| . 

TABLE II: Comparison of changes in multipole amplitudes 
resulting from addition of a sextupole element to the stan- 
dard fit model with Fourier coefficients of a SS ID Gaussian 
on azimuth. The comparison demonstrates equivalence be- 
tween the changes in multipole amplitudes (columns 2 and 3 
of Table |TJ) and the Fourier coefficients of a ID Gaussian on 
azimuth (column 5 of this Table). Column 5 is obtained by 
multiplying Fourier coefficients F m by the value Aid = 0.25 
that best matches the third column. That value corresponds 
to the difference AAi = 0.76 - 0.51 from Table U 



m 


parameter 


data 




F m A 1D 


1 


-AA D 


0.115 


0.445 


0.11 


2 


aa q 


0.05 


0.20 


0.05 


3 


AA S 


0.012 


0.053 


0.013 



Table |TT] (third column, data) shows the changes in 
fitted multipole amplitudes between the standard model 
and that including the sextupole term. The fourth col- 
umn, with F m defined in Eq. ([5]), shows the calculated 
Fourier coefficients for a unit-amplitude ID Gaussian on 
azimuth with width cf^^^d — 0.73 [16]. The fifth col- 
umn {F m A\D) shows the Fourier coefficients for a SS ID 
Gaussian with amplitude Am — 0.25. The Gaussian 
amplitude and width were adjusted to achieve the best 
match between columns 3 and 5, achieving equivalence 
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within small data uncertainties. The amplitude Aid of 
the added SS ID Gaussian on azimuth corresponds to 
the difference AAi between SS 2D peak amplitudes in 
Table U for the two fit models. And the inferred ID peak 
width approximates the azimuth width of the SS 2D peak 
in the standard fit model. Thus, the parameter changes 
resulting from inclusion of a sextupole element in the 2D 
model are equivalent to adding a constant offset Aid to 
the H(t]a) factor of the SS 2D peak model. 

Table U (fourth column) confirms the equivalence with 
a 2D fit that replaces the added sextupole element by 
an offset in the factor of the SS 2D peak model 
[A 1 exp(-r ? i/2 0; 2 A ) A' 1 exp(-r ] l/2a'l)+A 1D ]. The 
added constant in H(n^) then represents a SS ID Gaus- 
sian on (j>A in the 2D model [see Eq. (JTJ)]. The fit results 
are consistent with the previous exercise: the fitted ID 
Gaussian amplitude is Am — 0.29 compared to 0.25, and 
the modified 2D amplitude A[ = 0.47 is consistent with 
A 1D (i.e., A\ + Am = A 1 = 0.76). Note that in column 4 
Ap and Aq have reverted to their values for the standard 
fit model, since the added SS ID Gaussian on azimuth is 
modeled explicitly by a Gaussian function rather than by 
a truncated Fourier series. 




FIG. 8: Left: The curves represent the changes in three 
Fourier amplitudes (multipoles) between the standard fit 
model and a model including an added sextupole element 
applied to the same 2D data histogram, as reported in Ta- 
ble UU Subscript X represents dipole D, quadrupole Q and 
sextupole S for m = 1, 2 and 3 respectively. Right: The sum 
of the three sinusoids in the left panel per the axis label (solid 
curve) compared to a periodic ID Gaussian (dashed curve). 



Figure [5] illustrates the equivalence graphically. The 
left panel shows the changes in three multipoles, be- 
tween the standard 2D model and that with the added 
sextupole element. The sinusoid amplitudes are as in Ta- 
ble |TT] (third column). The right panel shows (solid curve) 
the sum of the three sinusoids in the left panel. The 
dashed curve is a SS ID Gaussian with amplitude 0.255 
and width 0.73 demonstrating the accurate equivalence. 
Remaining small differences are consistent with the miss- 
ing octupole term (m = 4) in the truncated Fourier series. 



C. Model ambiguities and parameter significance 

The previous exercise demonstrated that the nonzero 
sextupole amplitude resulting from fits with 2D model 
(b) in Table U is actually associated with the jet-related 
SS 2D peak. We now return to the significance of 2D 
model-fit results. As noted in Sec. IV Bl the standard fit 
model of the SS 2D peak in Eq. (TTJ is effectively fac- 
torized to ID functions on and 0a- Whatever the 
actual peak structure H(j]a) we do observe that the nar- 
row ID Gaussian G(0a) is independent of 77A, and fac- 
torization accurately describes the SS 2D peak. We can 
therefore simplify the model comparison in Sec. IVII Al 
to a single ID Gaussian on t^a with parameters (A, a^) 
vs a ID Gaussian plus constant offset with parameters 
(A',cr', v A 1D ). 




-0.5 
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FIG. 9: Left: The data histogram is the 2D histogram in 
Fig. [4] (left panel) projected onto ID tja- The dashed curve 
through the data histogram is a single ID Gaussian with pa- 
rameters (A,a v )- The dash-dotted curve is the sum of a ID 
Gaussian (A',<r^) plus a constant offset A\d represented by 
the dotted line. The two curves cannot be resolved within the 
n acceptance limits. Right: The second difference of the data 
histogram in the left panel corresponds to its local curvature. 
The calculated curvature A/a^ of the ID Gaussian in the left 
panel (dashed curve) is the dash-dotted line in this panel. 

Figure [5] (left panel) repeats the projected ID his- 
togram in Fig. [4] (right panel) . The projected histogram 
is accurately modeled by a ID Gaussian (dashed curve) 
with amplitude A = 0.20 and width a v = 2.7. The ID 
peak symmetric about t?a = can also be described in a 
model-independent way by a Taylor series in even powers 
of ?7a. Two terms are systematically significant within 
the X] acceptance: a constant term (peak amplitude) and 
a quadratic term (peak curvature at the mode). The 
quartic term is not statistically significant. 

In the right panel the Gaussian curvature value 
A/a^ = —0.027 (dash-dotted line) is plotted together 
with the bin- wise second difference of the data histogram, 
comparing local statistical noise to the curvature signal. 
The fourth difference (quartic) would be overwhelmed 
by statistical fluctuations. The Gaussian curvature is 
consistent with the Taylor series. The number of model 
parameters thus matches the data structure (significant 
series coefficients) and the data are accurately described. 
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If a constant offset is added to the ID fit model the 
data can constrain only certain combinations of the pa- 
rameters. The constraints are A = A' + Am for the 
common peak amplitude and A/o? = A 1 ja'^ for the 
common peak curvature at the mode. The parameters 
(A' , cr^, Aid) are otherwise free to vary within those con- 
straints. If we choose the value A%d = 0.077 (dotted line 
in the left panel, equivalent to Aid = 0.29 for the unpro- 
jected 2D peak model) the data constraints determine 
that A' = 0.123 and a' rj = 2.12. A dash-dotted curve 
representing the sum of the modified ID Gaussian and 
the constant offset is superposed in the left panel upon 
the dashed Gaussian curve. The curves are not distin- 
guishable within the rj acceptance limits. 

The same exercise could be carried out with a range 
of Any values including zero (the standard data model). 
Neither the offset parameter Am nor the mathemati- 
cally equivalent sextupole amplitude As is systematically 
significant for data within the STAR TPC acceptance. 
Changes in fitted As, Ad and Aq are strongly correlated 
and must be combined to determine the actual changes 
in the 2D model. In the present case the effective changes 
in the model are negligible, even though some individual 
parameters may vary substantially. Considering individ- 
ual model parameters in isolation (e.g., ignoring strong 
covariances) can be very misleading. 

VIII. NONJET "HIGHER HARMONICS" 

The (jet-related) SS 2D peak must contribute to higher 
multipoles inferred from ID Fourier analysis. And we 
demonstrated that the SS peak can also contribute an 
apparent jet-related sextupole for 2D model fits that are 
insufficiently constrained by data. But can it be demon- 
strated that there is a significant contribution to higher 
multipoles not associated with the SS 2D peak? In other 
words, are nonjet higher multipoles necessary to describe 
measured 2D angular correlations? We must establish 
what nonjet multipole elements are necessary to describe 
2D angular correlation data, and what elements are not 
required by the data or are excluded. 



A. Necessity of the nonjet quadrupole element 

Figure [10] provides a demonstration that the nonjet 
quadrupole is a necessary element of the standard 2D 
fit model. The p t -integral data are from 9-18%-central 
Au-Au collisions. The standard fit model provides an ex- 
cellent description of those data, and a nonjet quadrupole 
amplitude Aq = 0.18 ± 0.01 is inferred for that central- 
ity [J] . The fit residuals for the standard fit model shown 
in panel (d) are consistent with statistical uncertainties. 
The other panels reveal the consequences of omitting the 
nonjet quadrupole element from the fit model. 

The best-fit model in panel (a) does not include a non- 
jet quadrupole element. The corresponding fit residuals 




FIG. 10: Model fit to p t -integral data from 9-18%-central 
200 GeV Au-Au collisions demonstrating the necessity of a 
nonjet quadrupole model element, (a) fitted model without 
quadrupole element, (b) data (c) fit residuals without nonjet 
quadrupole element, (d) residuals with standard model fit. 

in panel (c) reveal at least 80% of the correct Aq ampli- 
tude near 0a = 7r. To minimize x 2 with the improper 
model the fitting routine has compensated by overesti- 
mating the SS 2D peak amplitude by about 30% (0.76 
— > 1.0), leading to the SS minimum in the residuals at 
4>A = 0. The fit residuals in panel (c) demonstrate that 
the nonjet quadrupole, a quadrupole component indepen- 
dent of the SS 2D peak, is a necessary model element. 



B. Adding a nonjet sextupole component to data 

Figure [IT] explores the same question for higher mul- 
tipoles: Can we establish the necessity of a nonjet sex- 
tupole from data analysis. We determine the sensitivity 
of the fitting process to a nonjet sextupole by deliberately 
adding such a term to the 2D data and refitting with the 
standard 2D model having no nonjet sextupole element. 

In Fig. [TTJ a sextupole component has been added to 
the input data with As{2D} = 0.034. The reason for 
that value is explained below. Panel (a) shows the stan- 
dard model fit, panel (b) shows the modified input data, 
and panel (c) shows the resulting fit residuals, with four- 
times increased sensitivity. The standard-model elements 
do minimize \ 2 by reducing the apparent sextupole com- 
ponent well below what was inserted. The compensa- 
tion process is essentially a reversal of the process de- 
scribed in Sec. lVII 51 Changes in the AS dipole (increase) 
and quadrupole (decrease) interact with a change in the 
SS 2D Gaussian (increase). Since Ad represents an AS 
dipole the changes in both dipole and quadrupole are 
subtracted from the SS peak amplitude increase, leav- 
ing a negative SS sextupole component to compensate 
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FIG. 11: Model fit to p t -integral data from 9-18%-central 200 
GeV Au-Au collisions to which a nonjet sextupole component 
has been added, (a) Standard fit model, (b) data with sex- 
tupole component added (c) fit residuals, (d) corresponding 
inferred SS 2D peak. 



A. Terminology: jet-related vs nonjet structure 

Proper description and interpretation of 2D angular 
correlations relies on maintaining distinctions among sep- 
arate correlation components and model elements. In 
Sec. HID I wc defined "nonjet" and "jet-related" structures 
in 2D angular correlations. The terminology is based on 
data systematics for p-p and more-peripheral Au-Au col- 
lisions where a jet interpretation for those structures is 
most likely Q. We then maintain that language consis- 
tently over the full centrality range of Au-Au collisions. 
Alternative interpretations of some parts of "jet-related" 
structure are possible in more-central Au-Au collisions. 

In describing multipole elements v m we distinguish be- 
tween (a) v m obtained from ID Fourier fits applied to all 
angular correlations projected onto ID 0a and (b) v m 
obtained from fits to angular correlation data with a 2D 
model that may include one or more multipole elements 
Ax- Multipoles from case (a) must mix jet-related and 
nonjet structure as discussed in Sec.lVIl Multipoles from 
case (b) may represent jet-related or nonjet contributions 
and are the principal subjects of this study. 



some fraction of that added to the data. The input value 
As{2D} — 0.034 chosen for the figure corresponds to 
reducing the usual Aq{2D} = 0.18 value to zero. If 
a larger sextupole amplitude were introduced the fitted 
quadrupole value would be driven negative to compen- 
sate. If the quadrupole amplitude is constrained to be 
positive definite a larger input A${2D} value would re- 
sult in significant sextupole structure in the residuals. 

Panel (d) shows the SS 2D peak inferred by subtract- 
ing fitted AS dipole and NJ quadrupole (zero amplitude) 
model elements from the modified data histogram. In- 
stead of A2D = 0.76 inferred from the unmodified data 
the SS amplitude is now 1.73, and the r.m.s. 77 width is 
4.3 rather than 2.7. Addition of a small sextupole com- 
ponent to the data causes large changes in inferred SS 2D 
peak structure equivalent to adding a narrow ID Gaus- 
sian with amplitude approximately 1. We conclude that 
attempts to isolate nonjet higher harmonics, multipoles 
not associated with the SS 2D peak, cannot succeed. 



IX. DISCUSSION 

Conjectured triangular flow and other flow-related 
multipoles offer the possibility to reinterpret the SS 2D 
peak in angular correlations as a flow phenomenon rather 
than a jet manifestation, seeming to buttress claims for 
formation of a flowing, dense partonic medium with small 
viscosity in high-energy A-A collisions. In this study 
we have identified several problems with recent higher- 
harmonic-flow claims. We confronted supporting argu- 
ments for higher harmonic flows with observed properties 
of 2D angular correlation data as reported in Ref. [4J. 



B. ID Fourier analysis vs 2D angular correlations 

The possibility to interpret most Au-Au angular cor- 
relation structure in terms of flows has motivated intro- 
duction of multiple sinusoids to correlation data model- 
ing. The prototypic example was "elliptic flow" analysis 
based on the "event-plane" V2{EP} and two-particle cu- 
mulant «2{2} methods [27j . Both are equivalent to fit- 
ting all angular correlations projected onto ID </>a with 
a single cos(20a) function [lj| [2i|. Sinusoid analysis 
was then extended to include "triangular flow" i<3{2} by 
adding a cos(30a) term to ID model fits [8j. Generaliza- 
tion to a multi-term truncated Fourier series has recently 
emerged [111, |33j]. In Sec. I VII we emphasized multipoles 
inferred by projecting all angular correlations, including 
the SS 2D peak, onto ID azimuth and some consequences 
of isolating individual SS peak Fourier components. 

Of course any ID distribution on periodic azimuth can 
be represented accurately by a discrete Fourier series — 
the representation is in some sense trivial. But a ID 
Fourier series cannot describe the strong 77A dependence 
of 2D angular correlations, cannot describe correlations 
in p-p and more-peripheral Au-Au collisions that most 
likely correspond to jet production. And a single ID se- 
ries cannot discriminate among several data components 
that may arise from distinct correlation mechanisms. 

That some data may be described in terms of sinu- 
soids (whatever the fit quality) does not imply that a 
flow interpretation is justified. The nonjet quadrupole is 
strongly correlated with initial-state A-A geometry, but 
its systematics appear to contradict a hydrodynamic in- 
terpretation [HHHl]. 

The terms of a ID Fourier series are orthogonal. But 
orthogonality of Fourier series terms does not imply that 
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a ID Fourier series is a preferable data model for 2D an- 
gular correlations. There are potentially three sources for 
any sinusoid in a ID Fourier series: (a) the SS 2D peak, 
(b) the AS ID peak, (c) nonjet sinusoids not associated 
with (a) or (b). The SS 2D peak can be reduced to sev- 
eral azimuth multipoles, but with possibly large curva- 
tures on ryA [Hj]. Multiple orthogonal Fourier terms can 
be strongly correlated by coupling to a single resolved 
correlation component. And a single Fourier term may 
be coupled to multiple correlation components and their 
causative mechanisms (e.g., "flow" and "nonflow"). 

The relation of Fourier series to SS and AS periodic 
peak arrays is described in Sec. IV A of Ref. [26|. Each 
peak has its own Fourier series. If the two series are 
combined each term must represent at least two contribu- 
tions. For instance, the dipole term of a single ID Fourier 
series representing all projected 2D angular correlations 
/l£>{2} = Ao{SS} — Ad{AS} is the difference between 
two large numbers, the dipole component of the SS 2D 
peak (positive) and the AS dipole (negative) represent- 
ing jet structure in at least some A- A collisions [4j. The 
small difference has been attributed to "rapidity-even di- 
rected flow" and "global momentum conservation" [34j . 
Jet structure is thereby obscured. The net quadrupole 
term represented by A Q {2} = A Q {SS} + A Q {2D} in- 
cludes the substantial quadrupole component of the SS 
2D peak as well as the nonjet component. The sextupolc 
term is derived from the SS 2D peak as AgjS'S'}. In con- 
trast, the standard fit model described in Sec. [TTT] accom- 
modates all correlation structure on 77A and 0a- Each 
model element represents a single observed correlation 
component, facilitating unique interpretations. 



C. Significance of multipoles in 2D data models 

The main problem addressed by this study is the sys- 
tematic significance of sextupole and higher multipole 
amplitudes that may be inferred by adding cos(?710a) 
terms with m > 2 to the standard 2D angular correla- 
tion model. If one or more higher multipole elements are 
added to the standard 2D model are resulting nonzero 
amplitudes physically meaningful? In this study we have 
isolated two major issues: (a) additional model elements 
may exceed what is required by correlation structure 
measured within some detector 77 acceptance leading to 
large systematic uncertainties in certain model parame- 
ters, and (b) some fractions of several model multipoles 
may combine to form a peaked structure on azimuth that 
plays an unanticipated role in the 2D data model. 

In Sec. IVIll we demonstrated that the SS 2D peak for 
data obtained from more-central Au-Au collisions within 
a limited r\ acceptance may only constrain two model 
parameters on t]a- Addition of a sextupole amplitude to 
the standard fit model is equivalent to addition of a third 
parameter (an offset) to the SS peak model. The extrane- 
ous model parameter effectively amplifies relatively small 
and systematically insignificant variations in the data to 



appear as relatively large but still systematically insignif- 
icant variations in the extraneous parameter, but also in 
other parameters now coupled to it. Although individual 
model parameters may then seem unstable the overall 
model is not changed significantly. The problem is re- 
solved by removing the extraneous parameter (s) per the 
principle of parsimony or Ockham's razor. 

In Sec. IVIIIBl we demonstrated that whereas the non- 
jet quadrupole is a necessary and unique element of the 
standard fit model independent of the SS 2D peak a non- 
jet sextupole is not a necessary element. A NJ sextupole 
component cannot be isolated (with systematic signifi- 
cance) from the SS 2D peak structure. What is uniquely 
defined is the SS 2D peak relative to fitted AS dipole and 
NJ quadrupole within the standard fit model. The struc- 
ture of the SS peak is factorizable, and each factor can be 
represented by a functional form. That unique structure 
can then be compared to theoretical predictions in whole 
or in parts according to explicit hypothesis. 



D. Modeling data with extreme conditions imposed 

For almost all 2D angular correlations from RHIC col- 
lisions the data histograms within some limited ry accep- 
tance can be described by the "standard fit model" con- 
sisting of an away-side ID dipole, a nonjet quadrupole 
and a same-side 2D peak. We have shown that any 
"higher harmonics" must then originate with the SS 2D 
peak. But for a small fraction of the data, associated 
with what may be called "extreme conditions," angu- 
lar correlation structure may deviate substantially from 
the standard fit model. Such conditions include high-p t 
"trigger-associated" cuts, cuts favoring large p-p collision 
multiplicities, acceptance cuts emphasizing large 77 values 
and "ultra-central" A-A collisions (0-1% central). 

For p-p and more-peripheral Au-Au collisions the en- 
tire SS 2D peak is resolved within the STAR TPC an- 
gular acceptance and the correlation model is fully con- 
strained. The SS 2D peak model is unambiguous: a 2D 
Gaussian. The SS peak continues to be well resolved un- 
til well above the sharp transition in peak properties at 

^trans 

However, in more-central collisions the elon- 
gated SS 2D peak extends sufficiently far outside the 
TPC n acceptance that its modeling on 77A becomes am- 
biguous. When only a fraction of the SS peak is observed 
some peak properties (e.g., amplitude and curvature at 
the mode) are more reliably determined than others (e.g., 
higher moments, tail structure). If a greater fraction of 
the SS 2D peak becomes accessible within a larger detec- 
tor acceptance more peak properties can be determined, 
including possible non-Gaussian tails at larger 7/a- 

Substantial deviations of same-side and away-side peak 
structure from the standard fit model are observed at the 
RHIC and LHC for certain combinations of the extreme 
conditions listed above. The SS 2D peak acquires long 
tails on 77A, and the AS ID peak develops a substan- 
tial deformation. All such deviations can be described in 
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terms of changes in the SS 2D or AS ID peak model func- 
tions. In the former case the change appears restricted 
to the ID H(t]a) factor, and modified Gaussian or non- 
Gaussian peaked functions are available (e.g., Sec. IV-C 
of Ref. HI). 

But a trend has emerged to decompose the SS 2D peak 
into "jet- like" (small ija) and "ridge- like" (large tja) com- 
ponents and to model the latter with a ID Fourier se- 
ries in common with all other correlation structure (i.e., 
nonjet quadrupole and AS ID peak). Single Fourier 
terms may actually represent multiple physical mecha- 
nisms, including jet formation. However, the convention 
is adopted to interpret each term as a harmonic flow, 
including the m = 1 dipolc. In that system the connec- 
tion with established jet structure in peripheral collisions 
is broken, and there can be no consistent description of 
centrality evolution. Consistency could be restored by ac- 
knowledging that whatever the imposed conditions, mod- 
eling the SS 2D and AS ID peaks separately is essential. 
A single ID Fourier series representing two peaked struc- 
tures plus other nonjet components is not interpretable. 

X. SUMMARY 

The information contained in 2D angular correlation 
data from nuclear collisions, both p t -integral and pt- 
differential, is extensive and could be used to eliminate 
invalid theoretical models and suggest valid physical in- 
terpretations. However, some analysis methods minimize 
and distort that information by (a) projection onto sub- 
spaces and (b) many-to-one and one-to-many mappings 
between correlation components and model elements. 

Recently, a major initiative has emerged to infer 
"higher harmonic flows" (azimuth multipoles) from cor- 
relation data. Fourier analysis of angular correlations 
projected onto ID azimuth from limited pseudorapid- 
ity intervals is used to infer such multipoles. However, 
ID Fourier analysis conflicts with 2D model fits to un- 
projected angular correlations in which apparent jet-like 
structure appears to dominate correlations. To help re- 
solve the impasse we have conducted a detailed study of 
2D angular correlation analysis methods and interpreta- 
tions relating to the question of "higher harmonic flows." 

A ID Fourier analysis of 2D angular correlations pro- 
jected onto azimuth must describe the periodic data ac- 
curately by definition. But there is nothing to compel in- 
terpretation of the series terms as flows. The unprojected 
2D structure can be decomposed precisely into several 
components that retain their identity from p-p collisions 



to central Au-Au collisions. One of those components, 
the same-side 2D peak, is strongly structured on pseu- 
dorapidity. A single ID Fourier series cannot model that 
component and therefore fails as a data model. More- 
detailed examination establishes that multiple correla- 
tion components contribute to some individual Fourier 
terms (many-to-one) and some correlation components 
contribute to multiple Fourier terms (one-to-many) lead- 
ing to severe interpretation difficulties. Further analy- 
sis establishes that "higher harmonics" inferred from ID 
Fourier analysis are at least dominated by the SS 2D 
(jet-related) peak. 

A potential remedy for such ambiguities might be 
achieved by adding higher multipole terms to the "stan- 
dard fit model" of angular correlations that successfully 
describes p-p and more-peripheral Au-Au correlations. 
However, we demonstrate in this study that additional 
multipole elements can conspire to contribute to the 
jet-related SS 2D peak model rather than representing 
sought-after nonjet multipoles that might be interpreted 
as flows. For more-central Au-Au data adding m > 2 
terms to the standard data model leads to systemati- 
cally insignificant changes in the SS 2D peak model. Al- 
though some model parameters may change substantially 
the overall data model (sum of the elements) is quite sta- 
ble and provides a unique 2D data description. 

The last part of the study was determination of the ne- 
cessity for nonjet higher multipoles in the data descrip- 
tion. A known sextupole contribution was added to a 
data histogram and fitted with the standard 2D data 
model. The result was a change in the inferred SS 2D 
peak structure but no appearance of the sextupole com- 
ponent in the fit residuals. In contrast, omission of the 
nonjet quadrupole element from the fit model produced 
an equivalent component in the fit residuals. The nonjet 
quadrupole is thus a necessary and distinct model ele- 
ment. But we conclude that it is impossible to isolate 
higher nonjet multipoles from the SS 2D peak structure. 

It has been argued that at least in more-central A-A 
collisions some part or all of the SS 2D peak is a flow 
phenomenon. Its Fourier components might then be in- 
terpreted legitimately as flow manifestations. However, a 
general model-independent analysis program should con- 
front the intact structure of the SS 2D peak to determine 
the correct H(ji&) functional form vs systematic condi- 
tions. Theory tests must accommodate the SS 2D peak 
in its entirety, with all its related phenomenology, or fail. 
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